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ABSTRACT 

We study small-amplitude, nonlinear pulsations of uniformly and differentially rotat- 
ing neutron stars employing a two-dimensional evolution code for general-relativistic 
hydrodynamics. Using Fourier transforms at several points inside the star, both the 
eigenfrequencies and two-dimensional eigenfunctions of pulsations are extracted. The 
centrifugal forces and the degree of differential rotation have significant effects on 
the mode-eigenfunction. We find that near the mass-shedding limit, the pulsations are 
damped due to shocks forming at the surface of the star. This new damping mechanism 
may set a small saturation amplitude for modes that are unstable to the emission of 
gravitational waves. After correcting for the assumption of the Cowling approximation 
(used in our numerical code), we construct empirical relations that predict the range 
of gravitational-wave frequencies from quasi-periodic post-bounce oscillations in the 
core collapse of massive stars. We also find that the fundamental quasi-radial mode 
is split, at least in the Cowling approximation and mainly in differentially rotating 
stars, into two different sequences. 



1 INTRODUCTION 

When neutron stars are formed as a result of a violent event 
(a core collapse, an accretion-induced collapse, or a binary 
neutron star merger), they are expected to pulsate non- 
linearly before various dissipative effects, such as viscosity, 
magnetic braking and shock formation, settle them down 
to a near-equilibrium configuration. All these violent events 
could lead to high rotational rates for the final objects, due 
to angular momentum conservation. Moreover, the final ob- 
ject is expected to rotate differentially, since the initial distri- 
bution of angular momentum has to be redistributed after 
collapse or merger. Several axisymmetric and nonaxisym- 
metric pulsation modes are expected to be excited during 
neutron star formation and could be important sources of 
gravitational waves. In addition, some modes could become 
unstable to the emission of gravitational radiation. The suc- 
cessful detection and identification of pulsation modes re- 
quires a detailed understanding of the eigenfrequency and 
eigenfunction of each mode-sequence, especially since pul- 
sation modes are high-frequency gravitational-wave sources. 
In addition, it is necessary to have a good understanding of 
the various damping mechanisms that will operate. 

There are several fundamental questions that have to 
be resolved before we can realistically hope for detection 
of gravitational waves from neutron star pulsations: Can 
neutron stars attain high angular momentum at birth? How 
strongly differential is the initial rotational profile? What 



kind of modes are excited? What is the initial amplitude for 
each mode? What is the frequency for each mode? Which 
mechanisms dampen the pulsations and on what timescale? 
To date, none of the above questions has a definitive answer, 
but some partial understanding is emerging. 

Initial rotation rate: Naively, the collapse of a rotat- 
ing stellar core should lead to an extremely rapidly rotating 
neutron star, just by the argument of angular momentum 
conservation. This picture is, however, made more compli- 
cated due to the presence of magnetic fields, wind-mass loss, 
centrifugal hangup during collapse, fall-back accretion and 
pre- or post-collapse binary interactions. In recent work by 
Heger et al. (2003) who consider evolutionary sequences of 
realistic rotating pre-collapse cores, it is suggested that the 
initial rotational period of a neutron star can indeed be of 
the order of 1ms or less, if the magnetic field is neglected. 
However, if the star passes through a red-supergiant phase, 
then the dynamo model for magnetic braking of the angular 
velocity (Spruit, 2002), increases the initial period of proto- 
neutron stars to several milliseconds. Wind-mass loss during 
a Wolf-Rayet/helium-star phase (Heger & Woosley, 2003) is 
another mechanism for slowing down the core of an evolved 
star. Nevertheless, more rapidly rotating proto-neutron stars 
could be obtained through other formation scenarios, such 
as fall-back accretion (see Watts & Andersson, 2002, and ref- 
erences therein) , pre-supernova binary interactions (Pfahl et 
al. 2002; Ivanova & Podsiadlowski, 2003; Podsiadlowski et 
al. 2003) and post-supernova accretion from a binary com- 
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panion (Langer et al. 2003) . A new scenario for millisecond- 
pulsar birth in globular clusters is the merger of binary white 
dwarfs (Middleditch, 2003). 

Recently, Villain et al. (2003) have studied in detail 
the quasi-stationary evolution from core collapse to the for- 
mation of a neutron star and find that in most cases an 
initially hot proto-neutron star contracts and spins-up on 
a time-scale of seconds, until reaching high rotation rates 
with T/\W\ > 0.1 - 0.2, (where T is the kinetic energy 
and W is the gravitational binding energy of the star). If 
the initial core is rapidly rotating, centrifugal hangup can 
occur during collapse. The rotational evolution of a proto- 
neutron star can also be affected by the details of the high- 
density equation of state (EOS hereafter), e.g., the appear- 
ance of hyperons in the core can yield a very rapidly ro- 
tating, but metastable, compact star (Yuan k, Heyl, 2003). 
Axisymmetric core collapse simulations in general relativity 
(Dimmelmeier, Font & Miiller, 2001, 2002a, 2002b; Shibata, 
2003) and in Newtonian gravity but with realistic initial 
conditions (Ott et al. 2003; Kotake, Tamada & Sato, 2003; 
see also Fryer et al. 2002 and references therein) have also 
shown the formation of rapidly rotating neutron stars. 

Differential rotation: The core collapse studies by 
Dimmelmeier et al. (2001, 2002a, 2002b), and a quasi- 
equilibrium treatment by Liu & Lindblom (2000) in the case 
of accretion-induced collapse of white dwarfs, have shown 
that rapidly rotating proto-neutron stars are born with a 
modest degree of differential rotation. A recent detailed 
analysis by Villain et al. (2003) shows that the length-scale 
on which the angular velocity changes within the star is of 
the order of 7-10 km. Hypermassive neutron stars created 
in a binary merger also have moderate differential rotation, 
with a somewhat shorter length-scale than in core collapse 
(see e.g., Shibata & Uryu, 2000). A fundamental question 
about the rotational profile is how quickly it will become uni- 
form. Shapiro (2000) and Cook, Shapiro & Stephens (2003) 
have suggested that magnetic braking could drive the star 
to uniform rotation on a timescale of only seconds. A more 
recent study by Liu & Shapiro (2003) shows that a magnetic 
field will cause turbulent motions in a differentially rotating 
star, driving it to uniform rotation. It has also been sug- 
gested that a differentially rotating, viscous proto-neutron 
star will be brought to uniform rotation due to turbulent 
mixing on a much shorter timescale (Hegyi, 1977), however, 
more detailed computations, including realistic composition 
gradients, are required in order to obtain better estimates. 

Mode excitation: In the rotational core collapse sim- 
ulations by Dimmelmeier et al. (2001, 2002a, 2002b) the 
proto-neutron star quickly settles into quasi-equilibrium af- 
ter core bounce. Still, pulsations are excited and survive for 
several oscillation periods. The dominant pulsation modes 
that are expected to be excited in rotational core collapse 
are the quasi-radial (I = m = 0) mode and the quadrupole 
(I = 2, m = 0) mode. Due to rotational couplings of non- 
radial terms in its eigenfunction, the quasi-radial mode be- 
comes a strong emitter of gravitational waves in rapidly ro- 
tating stars (it could, in fact, become the dominant mode 
in which gravitational waves are emitted). The amplitude 
of the oscillations in the density is estimated as several 
percent at the center of the star. Newtonian simulations 
(e.g., Monchmeyer et al. 1991; Zwerger & Miiller, 1997) 
also reached similar conclusions (for a recent review of core 



collapse simulations see New, 2002). In the simulations of 
binary neutron star mergers by Shibata and Uryu (2000, 
2002), quasi-periodic oscillations are excited, through which 
strong gravitational waves are emitted. The frequency of 
these oscillations suggest that they could correspond to spe- 
cific non-axisymmetric normal modes of the star. In addition 
to the non-axisymmetric modes, the axisymmetric quasi- 
radial mode could also be excited. 

Unstable Modes: In both proto-neutron stars and hy- 
permassive neutron stars created in binary mergers, gravi- 
tational radiation can drive several modes unstable, through 
the CFS mechanism (for a review, see e.g. Friedman & Lock- 
itch, 2001), provided the star is rotating sufficiently rapidly. 
In relativistic stars, the I = m = 2 /-mode becomes unsta- 
ble when T/|W| > 0.07 for uniform rotation (Stergioulas & 
Friedman, 1998; Morsink, Stergioulas & Blattnig, 1999) and 
at somewhat larger T/\W\ for differentially rotating stars 
(Yoshida et al. 2002). The I = m = 2 r-mode can become 
unstable at considerably lower rotation rates (for a review 
see Kokkotas & Andersson, 2001). The suppression of the 
r-mode instability by the presence of hyperons in the core 
(Jones, 2001; Linblom & Owen, 2002) is not expected to op- 
erate efficiently in rapidly rotating stars, since the central 
density is probably too low to allow for hyperon formation. 
Moreover, van Dalen & Dieperink (2003) find the contribu- 
tion of hyperons to the bulk viscosity to be two orders of 
magnitude smaller than previously estimated. If accreting 
neutron stars in Low Mass X-Ray Binaries (LMXB, consid- 
ered to be the progenitors of millisecond pulsars) are shown 
to reach high masses of ~ 1.8M , then the EOS could be 
too stiff to allow for hyperons in the core (for recent observa- 
tions that support a high mass for some millisecond pulsars 
see Nice, Splaver & Stairs, 2003). Both the /-mode and the 
r-mode, when unstable, will grow on a timescale of several 
seconds. 

Frequency of excited modes: Even though in a typi- 
cal nonrotating neutron star the fundamental radial and 
quadrupole modes have frequencies larger than about 1.5 
kHz, in the rapidly rotating models created in core collapse 
simulations, the frequency of these modes is significantly re- 
duced to well below 1 kHz (Dimmelmeier et al. 2002b), which 
is within the sensitivity window of current gravitational- 
wave detectors, such as VIRGO, GEO600, TAMA and 
LIGO. The large reduction in frequency is due to the fact 
that differential rotation allows much lower central densities 
than uniform rotation. In the case of a binary merger, the 
simulations by Shibata & Uryu (2002) have shown quasi- 
periodic oscillations of a few kHz. The r-mode, if unsta- 
ble, would emit gravitational waves at a high frequency of 
4/3 the rotational frequency. These sources of gravitational 
waves are very interesting for the proposed wide-band dual 
sphere detector (Cerdonio et al. 2001). When the /-mode 
becomes unstable, it can have an initial frequency of several 
hundred Hz, which then reduces to zero as the star spins 
down and approaches the rotation rate at which the mode 
becomes stable again. Therefore, the /-mode would ide- 
ally sweep through the sensitivity window of current laser- 
interferometric detectors, provided it can grow to a suffi- 
ciently large nonlinear amplitude (Lai & Shapiro, 1995). 

Damping of pulsations: The stable pulsations excited in 
the relativistic and Newtonian rotational core collapse sim- 
ulations are seen to be strongly damped. In an isolated star, 
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quasi-radial modes are thought to be damped mainly by 
shear viscosity and gravitational radiation (on a timescale 
of at least thousands of oscillation periods). In principle, 
any pulsation of a magnetized star will also be damped by 
the magnetic field (see e.g. McDermott et al. 1984; Carroll 
et al. 1986). In the case of quasi-radial pulsations of proto- 
neutron stars, however, we do not have specific estimates for 
the corresponding damping time. Once excited, the quasi- 
radial mode would be an ideal long-term monochromatic 
source of gravitational waves. However, the strong damping 
seen in the above simulations reveals that a different mech- 
anism operates in a proto-neutron star, immediately after 
core bounce (Pons 2003). The strong damping is due to the 
presence of a high-entropy envelope, which surrounds the 
newly created neutron star immediately after its birth. The 
quasi-radial oscillations penetrate into the envelope and are 
not properly reflected by a sharp surface (as happens in the 
case of an isolated star). The envelope thus absorbs much 
of the initial pulsation energy on a dynamical timescale (a 
similar damping has been observed in nonlinear evolutions 
of unstable relativistic stars, forced to migrate to the sta- 
ble branch of equilibrium models, see Font et al. 2002). By 
the time the proto-neutron star cools down (on a timescale 
of several seconds) the initial pulsation amplitude has di- 
minished. Therefore, the expected gravitational-wave signal 
from pulsations in core collapse will be strongly damped. 

In order to model the expected signal, a detailed anal- 
ysis of the damping mechanism is required (taking into ac- 
count various factors, such as the EOS, the initial rotating 
stellar core etc.). The unstable /- and r-modes grow on a 
timescale of the order of several seconds or longer, by which 
time the initial pulsations have been dramatically damped. 
In the case of a binary merger, a high-entropy envelope will 
also form. Its extent will depend on the mass ratio of the 
two stars - an equal-mass irrotational binary could create 
only a small envelope, which will not dampen the pulsa- 
tions as quickly as in a core collapse, while an unequal-mass 
binary merger could result in a somewhat larger envelope 
being formed (see Shibata, Taniguchi & Uryu, 2003) and 
consequently in a stronger damping of quasi-periodic pulsa- 
tions. In any case, the quasi-periodic oscillations in a binary 
merger will be much more long-lived than in a core collapse. 

Pulsation modes of rapidly rotating stars, assuming uni- 
form rotation, have been computed by Yoshida and Eriguchi 
(1999, 2001), Font, Stergioulas & Kokkotas (2000), Font et 
al. (2001) and Stergioulas & Font (2001) in the Cowling 
approximation (in which spacetime perturbations are ne- 
glected). In full general relativity, Stergioulas & Friedman 
(1998) and Morsink, Stergioulas & Blattnig (1999) com- 
puted the onset of the I — m — 2 /-mode instability, 
while Font et al. (2002) computed the frequencies of the 
two lowest-order quasi-radial modes. For differentially ro- 
tating stars, the only existing computation is by Yoshida et 
al. (2002), who obtained the I — m — 2 /-mode frequencies 
for a moderate strength of differential rotation. Nonlinear 
effects of radial pulsations in nonrotating relativistic stars 
have been studied by Sperhake, Papadopoulos & Andersson 
(2001) (see also Sperhake, 2002). 

In the present paper, we study in detail two different 
sequences of uniformly and differentially rotating relativis- 
tic polytropes, and obtain the eigenfrequencies and eigen- 
functions of several pulsations modes in the Cowling ap- 



proximation. For a sequence of fixed central energy density 
the mode-frequencies are only weakly affected by the rota- 
tion rate. However, for a sequence of fixed rest mass, the 
mode frequencies continuously decrease as the rotation rate 
increases (and the central density decreases). An interest- 
ing result of our study is that the fundamental quasi-radial 
mode splits into two different modes, an effect which is more 
prominently seen in differentially rotating models. Based on 
the frequency spectrum of nonrotating and uniformly rotat- 
ing stars, this split is not expected. If it turns out not to 
be really a new mode, it could be due to the fact that in 
the Cowling approximation the energy and momentum con- 
servation are violated. For example, this violation leads to 
the appearance of an unphysical "fundamental" dipole mode 
in simulations of nonrotating stars (Font et al. 2001). An- 
other important outcome of our investigation is that when 
studying the excitation of quasi-radial pulsations in rapidly 
rotating models (which are near their mass-shedding limit) 
we find that the pulsations are damped because of mass- 
shedding in the equatorial region. This new damping mech- 
anism could set a severe limit on the saturation amplitude 
of unstable modes. 

The rest of the paper is organized as follows. In Sec- 
tion |21 we outline our computational method, while in Sec- 
tion |3] we describe the equilibrium properties of all initial 
models we subsequently evolve. In Section [I] the eigenfre- 
quencies and eigenfunctions of pulsating stars are presented 
while Section |5] focuses on the splitting of the fundamental 
mode. Sections |S| and [7| discuss the damping of pulsations 
due to mass-shedding and the implications of our results on 
the detection of gravitational waves, respectively. Finally, 
a summary and a discussion of our results is presented in 
Section H 



2 OUTLINE OF COMPUTATIONAL METHOD 

We study the axisymmetric pulsations of rapidly rotating 
relativistic stars by first constructing several sequences of 
uniformly and differentially rotating equilibrium models, as 
described in detail in Section [3] The equilibrium models as- 
sume a relativistic polytropic EOS of the form 

P = Kp^ N , (1) 
e = p + Np, (2) 

where p is pressure, e is energy density, p is rest-mass density, 
N is the polytropic index and K is the polytropic constant. 
Unless otherwise noted, we choose dimensionless units for 
all physical quantities by setting c = G — Mq — 1. 

The nonlinear time-evolutions of perturbed equilibrium 
models are carried out using the same numerical hydrody- 
namics code developed by Font et al. (2001, 2002). Suitable 
perturbations of selected equilibrium variables are added to 
the equilibrium model, in order to excite specific modes. In 
the absence of the true eigenfunction of a given mode, each 
perturbation is chosen so as to mimic the angular depen- 
dence of the eigenfunction of the corresponding mode of a 
slowly-rotating Newtonian star. Usually, this ensures that 
the chosen mode will dominate the time-evolution at least 
for the slower rotating models. However, since the perturba- 
tion is not exact, additional pulsation modes will be excited, 



4 Stergioulas, Apostolatos and Font 



especially for rapidly rotating models. For the / = modes 
a density perturbation is used, of the form 

Sp = ap c sin I 7T — I , (3) 



where p c is the central density of the star and r s (8) is the 
coordinate radius of the surface of the star. The constant a is 
the amplitude of the perturbation, which we normally take 
to be of the order of 1%. The I = 2 modes are excited, by 
perturbing the ^-component of the four- velocity as follows: 



ug = a sm 



.(*) 



sm V cos f 



(4) 



(see Font et al. 2001 for more details). 

The perturbed models are evolved in time with a 
two-dimensional general-relativistic hydrodynamics code, in 
which we have implemented a Godunov-type scheme based 
on the Marquina flux formula and the 3rd-order PPM re- 
construction (see Font et al. 2000, 2001; see also Font, 2003, 
for a review of Godunov-type schemes in general relativistic 
hydrodynamics). Keeping the spacetime fixed to the initial 
equilibrium state during the evolution corresponds to the 
Cowling approximation in perturbation theory. Below the 
mass-shedding limit (see Friedman, Ipser & Parker, 1986, 
for the precise definition of the mass-shedding limit in the 
case of rapidly rotating relativistic stars), we assume that 
the star remains isentropic by enforcing the EOS Q,©. 
Near the mass-shedding limit, when shocks form (see Sec- 
tion the adiabatic ideal fluid EOS is used instead 



p = (r - l)pe, 



(•») 



where F — 1 + 1 /N and e is the specific internal energy. No- 
tice that at the initial time the isentropic equilibrium models 
constructed with the polytropic EOS are consistent 

with the ideal fluid EOS 0. 

The time series of the evolved perturbations are Fourier 
analyzed, and the peaks in the corresponding spectra are 
identified with specific pulsation modes, starting from the 
nonrotating member of the sequence, where the pulsation 
frequencies are known from perturbation theory. As the ro- 
tation rate increases, it becomes increasingly more difficult 
to identify specific modes in the Fourier spectrum. For this 
reason, we also extract the eigenfunction for each peak in 
the Fourier spectrum and use it as an additional criterion to 
identify specific modes (see Section [IJ. 



3 EQUILIBRIUM MODELS 

Equilibrium models of rotating relativistic stars are con- 
structed using the numerical code rns (Stergioulas & Fried- 
man, 1995) which was extended to include differential rota- 
tion. Extensive tests of the accuracy of the code in the case 
of uniform rotation can be found in Nozawa et al. (1998) 
and in Stergioulas (2003) . The metric describing an axisym- 
metric relativistic star is assumed to have the usual form 



where the metric functions 7,p,w and a depend on the co- 
ordinates r and 9 only 1 . Axisymmetry enforces the specific 
angular momentum measured by the proper time of matter 
j = u u<f, to be a function of the angular velocity CI only. A 
different specific angular momentum 



3 — u <t> 



£ +P 



(7) 



is locally conserved during the phase of homologous collapse 
of a rotating star. The Rayleigh criterion for local dynamical 
stability to axisymmetric perturbations is 



dCl 



< 0. 



(8) 



The simplest common choice of the differential rotation law 
j — j(Cl) that satisfies the Rayleigh stability criterion (see 
Komatsu, Eriguchi & Hachisu, 1989a,b) leads to an angular 
velocity distribution of the form 



ci c -n = — 

A 2 



(CI- 



(Cl-u)) 2 r 2 



(9) 



where Cl c is the angular velocity on the rotational axis and 
A is a parameter with units of length, that determines the 
length scale over which the angular velocity varies inside the 
star (in the limit of A — * oo, uniform rotation is recovered). 

When constructing sequences of fixed rest mass, the ra- 
dius of the star can vary by more than a factor of two. If 
one would characterize the sequence by a fixed value of the 
parameter A, then different models would correspond to dif- 
ferent degrees of differential rotation. In order to ascribe the 
same degree of differential rotation to all models along a 
sequence, we follow Baumgarte, Shapiro & Shibata (2000) 
in normalizing the parameter A by the radius of the star. 
Thus, we define 



A= A/r e , 



(10) 



ds 



Kd(f>-ujdty+e 2 °'(dr+rde 2 ), (6) 



where A is now a dimensionless parameter and r e is the 
equatorial coordinate radius of the star. Notice that in the 
Newtonian limit, A is the radius of the cylinder, as a fraction 
of the radius of the star, where the angular velocity falls to 
one half of the central angular velocity. 

Our focus is on the effect of rotation on pulsation 
modes. Hence, we do not survey a broad range of suggested 
high-density equations of state, but rather choose a poly- 
tropic EOS with TV = 1 and K = 100, respectively. This 
choice corresponds to models of neutron stars having mass 
and radius similar to those constructed with a realistic EOS 
of average stiffness. We focus attention on two different se- 
quences of differentially rotating models (sequences A and 
B) and their uniformly rotating counterparts (sequences AU 
and BU). The equilibrium properties of all models are dis- 
played in Table [T] All sequences terminate at the same non- 
rotating model (thus, models AO, AU0, B0 and BU0 all co- 
incide) . 

The differentially rotating sequence A and its corre- 
sponding uniformly rotating sequence AU are characterized 
by a fixed rest mass Mo = 1.506Mq. Along sequence A, 
the degree of differential rotation is held fixed at A = 1. 



1 The metric function p in Eq. J§J should not be confused with 
the rest-mass density in Eq. IT1 . 
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Table 1. Properties of the four sequences of equilibrium models (A is a sequence of fixed rest mass Mo = 
1.506Mq with A = 1, AU is the corresponding sequence of uniformly rotating models, B is a sequence of 
fixed central rest mass density p c = 1.28 X 10~ 3 with A = 1, and BU is the corresponding sequence of 
uniformly rotating models). All models are relativistic polytropes with N = 1 and K = 100. The definitions 
of the various quantities arc given in the main text. Notice that all quantities are in dimensionless units with 
c = G = Mq = 1. 



model 


Sc 

(xl0~ 3 ) 


M 


R 




r p /r e 


n c 

(xl0~ 2 ) 


n e 

(xl0~ 2 ) 


T/\W\ 


AO 


1.444 


1.400 


9.59 


8.13 


1.0 


0.0 


0.0 


0.0 


Al 


1.300 


1.405 


10.01 


8.54 


0.930 


0.202 


0.076 


0.018 


A2 


1.187 


1.408 


10.40 


8.92 


0.875 


0.258 


0.098 


0.033 


A3 


1.074 


1.410 


10.84 


9.35 


0.820 


0.294 


0.113 


0.049 


A4 


0.961 


1.413 


11.37 


9.87 


0.762 


0.319 


0.123 


0.066 


A5 


0.848 


1.418 


12.01 


10.49 


0.703 


0.334 


0.130 


0.086 


A6 


0.735 


1.422 


12.78 


11.25 


0.643 


0.338 


0.134 


0.107 


A7 


0.622 


1.427 


13.75 


12.21 


0.579 


0.334 


0.134 


0.131 


A8 


0.509 


1.433 


15.01 


13.45 


0.513 


0.320 


0.130 


0.158 


A9 


0.396 


1.439 


16.70 


15.13 


0.444 


0.295 


0.122 


0.189 


A10 


0.283 


1.447 


19.03 


17.44 


0.370 


0.260 


0.110 


0.223 


All 


0.170 


1.456 


21.92 


20.30 


0.294 


0.218 


0.094 


0.260 


AUO 


1.444 


1.400 


9.59 


8.13 


1.0 


0.0 


0.0 


0.0 


AU1 


1.300 


1.404 


10.19 


8.71 


0.919 


1.293 


1.293 


0.020 


AU2 


1.187 


1.407 


10.79 


9.30 


0.852 


1.656 


1.656 


0.037 


AU3 


1.074 


1.411 


11.56 


10.06 


0.780 


1.888 


1.888 


0.055 


AU4 


0.961 


1.415 


12.65 


11.14 


0.698 


2.029 


2.029 


0.076 


AU5 


0.863 


1.420 


14.94 


13.43 


0.575 


2.084 


2.084 


0.095 


BO 


1.444 


1.400 


9.59 


8.13 


1.0 


0.0 


0.0 


0.0 


Bl 


1.444 


1.437 


9.75 


8.24 


0.950 


0.180 


0.067 


0.013 


B2 


1.444 


1.478 


9.92 


8.36 


0.900 


0.257 


0.094 


0.026 


B3 


1.444 


1.525 


10.11 


8.49 


0.849 


0.319 


0.116 


0.040 


B4 


1.444 


1.578 


10.31 


8.63 


0.800 


0.373 


0.134 


0.055 


B5 


1.444 


1.640 


10.53 


8.77 


0.750 


0.423 


0.150 


0.071 


B6 


1.444 


1.713 


10.76 


8.91 


0.700 


0.471 


0.165 


0.087 


B7 


1.444 


1.798 


11.01 


9.05 


0.650 


0.519 


0.179 


0.105 


B8 


1.444 


1.899 


11.26 


9.17 


0.600 


0.568 


0.192 


0.124 


B9 


1.444 


2.020 


11.50 


9.26 


0.550 


0.623 


0.205 


0.144 


BIO 


1.444 


2.167 


11.71 


9.27 


0.500 


0.689 


0.219 


0.165 


Bll 


1.444 


2.341 


11.80 


9.13 


0.450 


0.777 


0.236 


0.187 


B12 


1.444 


2.532 


11.64 


8.72 


0.400 


0.912 


0.258 


0.207 


BUO 


1.444 


1.400 


9.59 


8.13 


1.00 


0.0 


0.0 


0.0 


BUI 


1.444 


1.432 


9.83 


8.33 


0.95 


1.075 


1.075 


0.012 


BU2 


1.444 


1.466 


10.11 


8.58 


0.90 


1.509 


1.509 


0.024 


BU3 


1.444 


1.503 


10.42 


8.82 


0.85 


1.829 


1.829 


0.037 


BU4 


1.444 


1.543 


10.78 


9.13 


0.80 


2.084 


2.084 


0.050 


BU5 


1.444 


1.585 


11.20 


9.50 


0.75 


2.290 


2.290 


0.062 


BU6 


1.444 


1.627 


11.69 


9.95 


0.70 


2.452 


2.452 


0.074 


BU7 


1.444 


1.666 


12.30 


10.51 


0.65 


2.569 


2.569 


0.084 


BU8 


1.444 


1.692 


13.07 


11.26 


0.60 


2.633 


2.633 


0.091 


BU9 


1.444 


1.695 


13.44 


11.63 


0.58 


2.642 


2.642 


0.092 
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Figure 1. Mass vs. central energy density plot of all models in se- 
quences A, AU, B and BU. Open circles correspond to uniformly 
rotating models, while the crosses correspond to differentially ro- 
tating models. The two nearly horizontal sequences arc A and AU, 
while the two vertical sequences are B and BU. Also shown are 
the sequence of nonrotating models (solid line) , and the sequence 
of models rotating at the mass-shedding limit for uniform rota- 
tion (dashed line). Differential rotation allows equilibrium models 
well beyond the region allowed for uniformly rotating models. 




1.0 



Figure 2. Density stratification in the fastest differentially ro- 
tating model of the fixed rest-mass sequence A. The maximum 
density appears off-center, at x/r e = 0.32. In comparison, the 
shape of the nonrotating star of same rest mass is shown (dashed 
line), scaled by the equatorial radius of the rotating model. 



The values of Mo and A are chosen in order to represent 
a newly-born, differentially rotating neutron star. The an- 
gular velocity at the equator is roughly 1/3 to 1/2 of the 
central angular velocity, which is similar to the degree of 
differential rotation obtained in typical core collapse simu- 
lations (see Villain et al. 2003). The fastest rotating model 
in sequence A has a polar to equatorial coordinate axes ra- 
tio of only r p /r e = 0.294, a ratio T/\W\ = 0.26 and rotates 
close to, but still below, the mass-shedding limit. The cen- 
tral density is nearly an order of magnitude smaller than the 
corresponding nonrotating model, while the circumferential 
radius is more than twice as large. The uniformly rotat- 
ing sequence AU only reaches an axes ratio of 0.575, a ratio 
T/\W\ — 0.095, half the central density, and a 50% larger ra- 
dius than the corresponding nonrotating model. Model AU5 
is at the mass-shedding limit. 

On the other hand, the differentially rotating sequence 
B and its corresponding uniformly rotating sequence BU 
are characterized by a fixed central density p c = 1.28 x 10~ 3 
and fixed A = 1. Its fastest rotating member has rest mass 
of Alo = 2.79Mq and a gravitational mass of 2.53Mq, cor- 
responding roughly to a hypermassive neutron star created 
in a binary neutron star merger. The degree of differential 
rotation is in rough agreement with simulations by Shibata 
& Uryu (2000, 2002). The axes ratio for the fastest rotating 
model is 0.4 (the model being still below, but close to, the 
mass-shedding limit). Since all models in the sequence are 
compact, the radius R only increases by 21%. The corre- 
sponding uniformly rotating sequence only reaches an axes 
ratio of 0.58 at the mass-shedding limit (BU9) with an in- 
crease in radius R by 40%. Thus, we see that when con- 



sidering a sequence of fixed central density, the uniformly 
rotating models attain a larger equatorial radius than dif- 
ferentially rotating models, which tend to expand out of the 
equatorial plane, becoming torus-like. 

Fig. shows all constructed equilibrium models in a 
mass vs. central energy density plot. The two nearly horizon- 
tal sequences are A and AU, while the two vertical sequences 
are B and BU. Open circles correspond to uniformly rotat- 
ing models, while the crosses correspond to differentially ro- 
tating models. Also shown are the sequence of nonrotating 
models (solid line) and the sequence of models rotating at 
the mass-shedding limit for uniform rotation (dashed line). 
Differential rotation allows equilibrium models well beyond 
the region allowed for uniformly rotating models (see also 
Baumgarte, Shapiro & Shibata, 2000, where such configura- 
tions are called hypermassive) . 

In Fig. [5] the density stratification for the fastest rotat- 
ing model of sequence A (All) is shown. In the figure, Carte- 
sian coordinates, x = rsm6,z — rcos8 are used, scaled by 
the equatorial coordinate radius r e . The maximum density 
is attained off-center. The dashed line shows the spherical 
surface of the nonrotating model with the same rest mass 
(AO), scaled by the equatorial coordinate radius of the ro- 
tating model All. It becomes evident that the large rotation 
rate of model All causes the equatorial radius to increase 
by more than a factor of two. 

Finally, we note that apart from the above two se- 
quences, we have also investigated the effect of the degree 
of differential rotation on the axisymmetric pulsations, by 
constructing various sequences that differ only in the value 

of A. 
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Table 2. Frequencies of the fundamental quasi-radial (I = 0) 
mode, F, its first overtone, Hi, the fundamental quadrupole 
(I = 2) mode, 2 / and its first overtone, 2 pi, for the sequence 
of uniformly rotating models AU. 



model 


F 


Hi 


2 / 


V 




kHz 


kHz 


kHz 


kHz 


AUO 


2.706 


4.547 


1.846 


4.100 


AU1 


2.526 


4.246 


1.800 


3.862 


AU2 


2.403 


4.090 


1.744 


3.592 


AU3 


2.277 


3.937 


1.663 


3.265 


AU4 


2.141 


3.795 


1.547 


2.847 


AU5 


1.960 


3.647 


1.330 


2.560 



Table 4. Same as Table |2] but for the sequence of uniformly 
rotating models BU. 



model 


F 


Hi 


2 / 


V 




kHz 


kHz 


kHz 


kHz 


BU0 


2.706 


4.547 


1.846 


4.100 


BUI 


2.657 


4.467 


1.855 


4.040 


BU2 


2.619 


4.409 


1.860 


3.944 


BU3 


2.579 


4.385 


1.857 


3.814 


BU4 


2.535 


4.371 


1.844 


3.645 


BU5 


2.495 


4.356 


1.815 


3.456 


BU6 


2.456 


4.357 


1.762 


3.244 


BU7 


2.417 


4.337 


1.686 


3.010 


BU8 


2.328 


4.300 


1.588 


2.710 


BU9 


2.313 


4.280 


1.558 


2.642 



Table 3. Same as Table El but for the sequence of differentially 
rotating models A. Also shown are the frequencies of an addi- 
tional fundamental mode Fu that appears mainly in differentially 
rotating models (at least in the Cowling approximation). 



model 


F u 


F 


Hi 


2 / 


2 Pi 




kHz 


kHz 


kHz 


kHz 


kHz 


AO 


2.706 


2.706 


4.547 


1.846 


4.100 


Al 


2.485 


2.561 


4.310 


1.822 


3.961 


A2 


2.361 


2.480 


4.163 


1.780 


3.822 


A3 


2.243 


2.386 


4.029 


1.738 


3.642 


A4 


2.142 


2.295 


3.900 


1.677 


3.442 


A5 


2.039 


2.201 


3.748 


1.598 


3.211 


A6 


1.921 


2.101 


3.563 


1.510 


2.973 


A7 


1.779 


1.982 


3.329 


1.403 


2.711 


A8 


1.609 


1.846 


3.120 


1.274 


2.494 


A9 


1.424 


1.667 


2.857 


1.132 


2.185 


A10 


1.223 


1.422 


2.503 


0.966 


1.905 


All 


1.044 


1.220 


2.174 


0.820 


1.658 



Table 5. Same as Table El but for the sequence of differentially 
rotating models B. 



model 


Fu 


F 


Hi 




2 / 


2 Pi 




kHz 


kHz 


kHz 


kHz 


kHz 


B0 


2.706 


2.706 


4.547 


1 


.846 


4.100 


Bl 


2.627 


2.658 


4.446 


1 


.880 


4.102 


B2 


2.561 


2.637 


4.421 


1 


.900 


4.090 


B3 


2.525 


2.632 


4.405 


1 


.913 


4.045 


Bl 


2.506 


2.632 


4.403 


1 


.924 


3.983 


B5 


2.487 


2.632 


4.422 


1 


.929 


3.907 


B6 


2.459 


2.633 


4.436 


1 


.925 


3.828 


B7 


2.423 


2.635 


4.447 


1 


.909 


3.721 


B8 


2.394 


2.646 


4.444 


1 


.890 


3.632 


B9 


2.360 


2.653 


4.413 


1 


.867 


3.567 


B10 


2.330 


2.662 


4.360 


1 


.842 


3.500 


Bll 


2.318 


2.678 


4.300 


1 


.832 


3.470 


B12 


2.328 


2.722 


4.240 


1 


.830 


3.473 



4 EIGENFREQUENCIES AND 
EIGENFUNCTIONS 

Although we use a nonlinear evolution code to study pul- 
sations of rotating stars, we restrict attention to small- 
amplitude pulsations (small in the sense that e.g. Sp/p ~ 
10 -2 ). Therefore, the time-evolution of a perturbed star can 
still be viewed (to a good approximation) as a superposition 
of linear normal modes. When obtaining the Fourier trans- 
form of the time-evolution of several variables, we verify that 
the various modes that are excited have indeed discrete fre- 
quencies (same frequency at any point inside the star in the 
coordinate frame). Thus, even though the evolution is non- 
linear, the amplitude is small enough to justify the use of 
the terms "eigenfrequency" and "eigenfunction" for the var- 
ious pulsation modes. In order to compute the real part of 
the eigenfrequency of a pulsation mode we Fourier-transform 
the time series of the evolution of a suitable physical vari- 
able (the density for the / = modes and vg for the I — 2 
modes). Instead of examining the Fourier spectra at a few 
specific points inside the star, we integrate the amplitude of 
the Fourier transform along a coordinate line, e.g. for the 
I = modes we examine the integrated Fourier amplitude 
along 6 = 7t/2 (equatorial plane), while for the I = 2 modes 



the integrated Fourier amplitude along a line of 9 = 7r/4 is 
used. 

Since the trial-eigenfunction used for exciting the pul- 
sations does not correspond exactly to a particular mode, 
several additional modes are excited, apart from the main 
mode one wishes to study. This is particularly true for very 
rapidly rotating models, where rotational coupling effects 
are significant and higher-order coupling terms in the mode- 
eigenfunctions become comparable to the dominant term. 
Thus, in order to identify specific modes for several models 
along a sequence, one has to begin with pulsations of a non- 
rotating star (for which the modes are identified by compar- 
ison to results obtained with linear perturbation theory, see 
Font et al. 2000, 2001, 2002) 2 and gradually identify modes 
for more rapidly rotating stars by comparing the peaks in 
a Fourier transform to the corresponding peaks in a more 
slowly rotating model. Close to the mass-shedding limit, 
avoided crossings between low-order and high-order modes 
can complicate the picture and lead to erroneous identifica- 
tions. For this reason, we do not only rely on the Fourier 
transforms at a few points inside the star, but reconstruct 

2 For a review of pulsations of nonrotating relativistic stars see 
Kokkotas & Schmidt (1999). 
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Figure 3. Amplitude of the Fourier transform of the time- 
evolution of the density in the equatorial plane, for model B7, 
after applying an I = perturbation. Several excited pulsation 
modes can be identified. For a particular mode, the amplitude of 
the Fourier transform correlates with its eigcnfunction. 



the whole two-dimensional eigenfunction of each mode, us- 
ing Fourier transforms at every point inside the star. At 
the eigenfrequency of a specific mode, the amplitude of the 
Fourier transform correlates with its eigenfunction. A change 
in sign in the eigenfunction corresponds to both the real 
and imaginary part of the Fourier transform going through 
zero. Comparing the eigenfunctions corresponding to differ- 
ent peaks in a Fourier transform, allows for an unambiguous 
identification of specific mode sequences. 

A example of the eigenfunction extraction using Fourier 
transforms is shown in Figure [3] which displays the ampli- 
tude of the Fourier transform of the time-evolution of the 
density, after an I — perturbation was applied to model 
B7. In the range of frequencies shown in Figure the am- 
plitude of the Fourier transform clearly correlates with the 
eigenfunctions of the fundamental F- and 2 /-modes (the ad- 
ditional mode, Fn, is discussed in Section |^J. Notice that in 
displaying the eigenfunctions we use the amplitude of the 
Fourier transform, multiplied by the sign of its real part. 

When comparing the eigenfrequencies between models 
that differ only in the degree of differential rotation, we see 
that a moderate degree of differential rotation has some (but 
not dramatic) effect on the eigenfrequencies. When A is de- 
creased significantly below ~ 1, which leads to a "strongly" 
differentially rotating model, most of the angular momen- 
tum of the star is concentrated in a narrow region around 
the rotational axis. Outside this region the star is only slowly 
rotating and the eigenfrequencies of its pulsations become 
again similar to those of a nonrotating model. Thus, the ef- 
fect of differential rotation on the eigenfrequencies becomes 
strongest, not for very small values of A (strong differential 



rotation), but for values that correspond to stars with only 
a moderate degree of differential rotation. 

4.1 Fixed rest mass sequences 

It is well-known that the frequencies of the fundamental 
I — and I — 2 polar modes of oscillation depend mainly 
on the central density of a star, or, equivalently, on the 
compactness M/R (see e.g. Hartle & Friedman 1975). This 
is particularly true for the axisymmetric (m = 0) modes. 
The sequences of fixed rest mass Mo = 1.506Mq start with 
a nonrotating model with compactness M/R = 0.15. Ro- 
tation increases the radius and decreases the central den- 
sity. The uniformly rotating sequence AU terminates at the 
mass-shedding limit, with a compactness of M/R = 0.095. 
The differentially rotating sequence can reach higher rota- 
tion rates and terminates near the mass-shedding limit with 
a compactness of M/R — 0.066. Based on the significant 
decrease of the compactness along the fixed-rest-mass se- 
quences, we expect a corresponding decrease in the frequen- 
cies of the fundamental modes (and a similar tendency for 
the first overtones). 

Table [5] displays the computed eigenfrequencies for the 
fundamental and first overtone of the I — and I = 2 
modes for the sequence AU. Rotation reduces the frequency 
of the fundamental quasi-radial mode from 2.71 kHz to 
1.96 kHz with corresponding changes in the frequencies of 
the other modes. Figure 2] displays the variation of the 
mode- frequencies with increasing T/|W| along sequence AU 
(dashed lines) . The rate of decrease in frequency for the fun- 
damental quasi-radial modes, with increasing rotation rate 
becomes larger as the mass-shedding limit is approached 
(due to a cusp forming in the equatorial region), while the 
decrease in frequency for the other modes is a nearly linear 
function of T/\W\. 

Table |3 reports the corresponding eigenfrequencies for 
the differentially rotating sequence A, which are also shown 
in Figure [I] (solid lines) . The frequencies of all modes de- 
crease nearly linearly with increasing T/|W|. Due to the dif- 
ferential rotation, the outer layers of the star rotate slower 
and the equatorial radius is smaller than a uniformly ro- 
tating model of same T"/|W|. This leads to a smaller sound- 
crossing time and correspondingly higher fundamental mode 
frequencies for the differentially rotating models. This ex- 
plains why the lines in Fig. 0] corresponding to the funda- 
mental modes of sequence A have smaller slopes than those 
corresponding to sequence AU. 

For the fastest rotating model of sequence A, the fun- 
damental quasi-radial mode has a frequency of only 1.22 
kHz, with the fundamental quadrupole mode having a fre- 
quency of 0.82 kHz. It should be emphasized that the above 
frequencies are computed in the Cowling approximation, 
which leads to higher values than the actual frequencies 
(see Yoshida & Kojima, 1997). Therefore, the actual fre- 
quencies of the fundamental quasi-radial modes should be 
in the range of 40 — 45% smaller than those computed here. 
The actual frequencies of the fundamental quadrupole mode 
should be roughly 15 — 20% smaller than our values. For the 
nonrotating model of sequences A and AU, the actual fre- 
quencies of the fundamental I — and I — 2 modes are 1.44 
kHz and 1.58 kHz, respectively (Font et al. 2002). Thus, 
the fundamental I = mode has a frequency only some- 
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Figure 4. Eigenfroquencies of various modes for sequences A and 
AU. The lower of the two lines that start at the F-mode frequency 
for the nonrotating models is an additional mode, Fn, appearing 
mainly in differentially rotating stars (at least in the Cowling 
approximation) . 



what smaller than the fundamental 1 = 2 mode (while in 
the Cowling approximation it has a frequency significantly 
larger than the fundamental I = 2 mode). 

In Font et al. (2002) it was found that, in the case of 
uniform rotation, both the actual frequency and the fre- 
quency in the Cowling approximation of the fundamental 
quasi-radial mode decrease in a similar way as the rotation 
rate increases. This leads to the following empirical rela- 
tion between the actual frequency and the frequency in the 
Cowling approximation. 



fn = f ( n c) + ( f - n 



(C) 



(11) 



where the superscript (C) denotes the Cowling approxima- 
tion, /n is the frequency of the fundamental quasi-radial 
mode for a rotating star with angular velocity Q, and fo 
is the corresponding frequency in a nonrotating star. This 
empirical relation was found to be accurate to within 2% at 
all rotation rates, even near the mass-shedding limit. Since 
the frequencies of differentially rotating stars are not much 
different with respect to those of uniformly rotating stars 
with the same oblateness, the above relation should approx- 
imately hold for differentially-rotating models as well. 

Based on Yoshida & Kojima (1997) and Yoshida & 
Eriguchi (2001) we can estimate the error in the fundamen- 
tal mode frequencies, due to the Cowling approximation. For 
the fastest rotating model of sequence A, with M/R ~ 0.07, 
we estimate that the actual fundamental I = and I = 2 
frequencies are 63% and 79% of our results in the Cowling 
approximation, respectively (actual frequency meaning the 
frequency without the assumption of the Cowling approx- 
imation). Taking into account the nearly linear scaling of 
the frequencies with increasing rotation rate in Figure^] we 
construct the following empirical relation 
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Figure 5. Same as Figure El but for the sequences B and BU. 



/(kHz) 



1.44 - 2.59- 



W\ 



(12) 



for the actual frequency of the fundamental quasi-radial 
mode and 



/(kHz) 



1.58 - 3.69- 



\W\ 



(13) 



for the actual frequency of the fundamental quadrupole 
mode for the differentially rotating models of sequence A. 
We estimate the uncertainty of the above empirical relations 
to be on the order of a few percent. 

According to the above empirical relations, for rapidly 
rotating proto-neutron stars with T/|W| in the range 0.14 — 
0.26, the frequency of the two fundamental modes will be 
in a range of roughly 0.65 — 1.1 kHz. This range of frequen- 
cies agrees well with the frequencies of gravitational waves 
observed in the rotating core collapse simulations by e.g. 
Dimmelmeier et al. (2001, 2002a, 2002b), which confirms 
the validity of our chosen equilibrium models and pulsation 
modes to model the gravitational waves produced in the 
above simulations 3 . 

The additional fundamental frequency Fu in Tables [5] 
and Eland in Figure 2] is discussed in Section 



4.2 Fixed central density sequence 

Another sequence of models with differential rotation that 
has been analyzed is sequence B, consisting of twelve models 
with the same central energy density e c = 1.28 x 10~ 3 but 
with axes ratios ranging from 1.0 to 0.40. The most rapidly 
rotating model corresponds to a hypermassive neutron star, 
such as those created temporarily in a binary neutron star 



3 Here we only refer to the regular collapse cases in Dimmelmeier 
et al. and not to multiple-bounce cases. Notice that, even though 
our polytropic index is different than in Dimmelmeier et al., our 
choice of the polytropic constant leads to models of similar com- 
pactness, for the same mass. 
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Figure 6. Two-dimensional eigenfunction of the density pertur- 
bation, corresponding to the I = 2 fundamental mode for the 
nonrotating model BO. The eigenfunction is shown in Cartesian 
coordinates with x and z being scaled by the equatorial coordi- 
nate radius r e , while 8p is in arbitrary units. 



merger. The corresponding sequence of uniformly rotating 
models, BU, of same fixed central energy density, comprises 
models with axes ratios from 1.0 to 0.58. The computed 
frequencies of the / = 0, 2 fundamental modes and their first 
overtones are shown in Figure and in Tables 2] and (the 
additional fundamental mode Fu is discussed in Section |SJ. 

The frequencies of the fundamental 1 = and 1 = 2 
modes changes by as much as 16% along the sequence BU. 
In comparison, for sequence B the eigenfrequencies for both 
the / = and / = 2 fundamental modes remain relatively 
unaffected by the rotation rate, changing by less than 3% 
and 5%, respectively, compared to their values in the non- 
rotating limit. This leads to the conclusion that a moderate 
amount of differential rotation, of the order of A = 1, makes 
the fundamental frequencies relatively insensitive to rotation 
and depending mainly on the central energy density of the 
star. This finding should simplify attempts to extract infor- 
mation about the physical properties of neutron stars in bi- 
nary mergers, when such events become observable through 
their gravitational- wave emission. 

As discussed in Section|2]we obtain the two-dimensional 
eigenfunction of a specific mode by computing Fourier trans- 
forms of selected variables at every point inside the star. 
An example of such an eigenfunction is shown in Figure |S] 
which displays the eigenfunction corresponding to the den- 
sity perturbations due to the 1 = 2 fundamental mode, for 
the nonrotating model BO. Since the model is nonrotating, 
the eigenfunction shows the expected quadrupolar structure. 
For the most rapidly rotating model B12, however, the eigen- 
function for the same mode appears severely altered (Figure 
with respect to the nonrotating limit. There are two main 
effects caused by rapid rotation. Firstly, the I = 2 funda- 
mental mode does not only couple to higher order spherical 
harmonics, but also to an I = term. This causes a signifi- 
cant density variation at the center of the star, whereas for 




Figure 7. Same as Figure [S] but for the most rapidly rotating 
model of sequence B (B12). 

the nonrotating model the center of the star is not pulsating. 
Secondly, the centrifugal force weakens the effective gravity 
in the equatorial region, causing the eigenfunction to attain 
a large amplitude near the equatorial surface of the star. 
The effect of the centrifugal force becomes extreme when 
the star rotates at the mass-shedding limit, as discussed in 
Section |B] The above effects are common in other modes, 
too, and also in all four sequences considered. 



5 SPLITTING OF THE FUNDAMENTAL 
MODE 

An interesting feature visible in Figures 2] and |K| is the ap- 
pearance of an additional fundamental mode, which we de- 
note as Fu. This mode is degenerate in the nonrotating limit, 
but becomes distinct from the F-mode for rotating mod- 
els. For sequence A, the Fn-mode has a similar dependence 
on rotation rate as the F-mode, with the difference in fre- 
quency between the two modes remaining nearly constant 
with increasing rotation rate. For sequence B, however, the 
frequency of the Fn-mode is decreasing with increasing ro- 
tation rate, while the frequency of the F-mode is increasing. 

For both sequences A and B (which have a fixed 
strength of differential rotation, A = 1) the Fn-mode is ex- 
cited at the same level as the F-mode, by the perturbation 
shown in Eq. J3J (the amplitude of the Fourier transform is 
similar for both modes). Varying the degree of differential 
rotation (for the same model and the same applied pertur- 
bation) we find that, while the amplitude of the F-mode 
remains constant, the amplitude of the Fn-mode decreases 
as the star becomes less differentially rotating. In the limit 
of uniform rotation, there still exists a peak in the Fourier 
transform, corresponding to the Fu mode, although it has 
an amplitude much smaller than the main F-mode (see Fig- 
ure|H] which compares the Fourier transforms for models B7 
and BU7). Extracting the eigenfunction of the radial velocity 
perturbation, in the equatorial plane, for model BU7, shows 
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Figure 8. Comparison of integrated Fourier amplitudes of the 
density evolution for a differentially rotating (B7) and a uniformly 
rotating (BU7) model. The peaks corresponding to the fundamen- 
tal quasi-radial mode F, the additional fundamental mode Fn and 
the fundamental 1 = 2 mode are shown. 
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Figure 9. Comparison of the eigenfunctions, in the equatorial 
plane, corresponding to variations in the velocity component v r 
for the fundamental quasi-radial mode F and for the additional 
fundamental mode F\\ (for the uniformly rotating model BU7) 



that the eigenfunction of the Fn-mode has a node close to 
70% of the equatorial radius, while the F-mode clearly has 
no node (see Figure . Even though the Fn-mode has a 
node in the radial velocity, which is a property of a first 
overtone in a nonrotating star, we still refer to it as a fun- 
damental mode, due to its degeneracy with the -F-mode in 
a nonrotating star. 

The distinction between the F-mode and the Fh-mode 
becomes very clear in the extracted two-dimensional eigen- 
functions. Figures [TU1 and [Til show the eigenfunction of the 
two modes, corresponding to the density perturbation, for 
the differentially rotating model B8. The eigenfunction of 
the F-mode is significantly modified by rotational couplings 
near the equatorial surface, but otherwise it is similar to the 
expected eigenfunction of the F-mode in a nonrotating star. 
In contrast, the eigenfunction of the Fn mode is similar to 
an F-mode eigenfunction in a nonrotating star only in the 
central region, but otherwise it is very different in the polar 
and equatorial regions. The shape of the eigenfunction in the 
equatorial plane (having no node in the density perturba- 
tion) is reminiscent of density perturbations in differentially 
rotating tori. 

Is the Fn-mode physical? Since our numerical code im- 
plements the Cowling approximation we cannot answer this 
question at present. In the Cowling approximation the en- 
ergy and momentum conservation are violated. In previ- 
ous simulations (Font et al. 2001) this violation of the con- 
straints has lead to the appearance of an unphysical "fun- 
damental" dipole mode. It is possible that the Fn-mode ob- 
served in the present simulations is also an artifact of the 
Cowling approximation. This can only be confirmed by new 
simulations in full general relativity. 



6 MASS-SHEDDING-INDUCED DAMPING OF 
PULSATIONS 

Linear perturbations of rotating stars are assumed to have 
a vanishingly small amplitude, so that the background equi- 
librium star is unaffected by a linear oscillation mode. How- 
ever, when one considers a finite-amplitude oscillation, then 
nonlinear effects can become important. Our nonlinear evo- 
lution code allows us to investigate such effects. 

At low rotation rates, a small oscillation amplitude of 
the order of 10 -2 does not lead to significant nonlinear ef- 
fects. But, as the star approaches the mass-shedding limit, 
the effective gravity near the equatorial surface diminishes. 
Exactly at the mass-shedding limit fluid elements are only 
marginally bound to the surface of the star. A small radial 
pulsation then suffices to cause mass-shedding after each os- 
cillation period. At an oscillation frequency around 2.3 kHz 
(in the Cowling approximation), matter that has been shed 
from the star does not have the time to fall back before more 
matter is shed (this statement applies to the specific exam- 
ple with an initial pulsation amplitude of the order of 10 -2 , 
that is studied here - more generally, the fall-back time will 
depend on the initial amplitude). In this way, a low-density 
toroidal envelope is created in the equatorial region, which 
expands with every oscillation period. 

Figure [T2"l shows the profile of the specific internal en- 
ergy e at two different times, for the uniformly rotating 
model BU9, which is at the mass-shedding limit. It is clearly 
seen that the matter is shed in the form of high-entropy 
shock waves. A first shock wave, shown at t=0.14ms, leaves 
the star as a result of the initial perturbation applied to the 
equilibrium model. After one oscillation period, a second 
shock wave leaves the star and both shock waves are shown 
at t=0.92ms. In the same fashion, consecutive shock waves 
are created after each oscillation period and the density in 
the toroidal envelope gradually increases. 
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Figure 10. Two-dimensional eigenfunction of the density per- 
turbation, corresponding to the / = fundamental mode for the 
differentially rotating model B8. 
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Figure 11. Same as Figure lTol but for the additional fundamen- 
tal mode Fu (see discussion in Section IS! . 



Since matter is shed from the star in the form of shocks, 
they carry away kinetic energy, to the expense of the pul- 
sational energy. In this way, the pulsations of the star are 
gradually damped. Figure 1131 shows a comparison between 
the time-evolution of the central rest-mass density in the 
slowly rotating model BUI and model BU9. It is evident 
that the damping of the pulsations due to mass-shedding 
is very strong. The damping will also be effective in stars 
which rotate at somewhat smaller rotation rates. The pre- 
cise damping rate will depend on several factors, such as 
the amplitude of the pulsation, the rotation rate of the star, 
the EOS etc. The assumption of the Cowling approximation 
which we have adopted in our simulations will also affect 
the damping rate, as it causes an imbalance between the 
different forces acting on a fluid element on the equatorial 
surface of the star. Simulations which take into account the 
time-evolution of the gravitational field during the oscilla- 
tions should yield more accurate damping rates. 

Our finding of the mass-shedding induced damping of 
pulsations in critically rotating stars can have severe conse- 
quences for unstable pulsation modes, such as /-modes and 
r-modes, since, over a timescale of several seconds, even a 
very small damping rate could suffice to saturate their am- 
plitude at values much less than order unity. 



7 DETECTABILITY OF GRAVITATIONAL 
WAVES 

Pulsating rotating neutron stars are gravitational-wave 
sources that depend on several parameters (EOS, mass, an- 
gular momentum, differential rotation law, initial amplitude, 
damping mechanisms etc.). All these parameters may have 
different effects on the oscillation spectrum of the star and, 
therefore, the successful extraction of the physical character- 
istics of the source from the gravitational-wave signal will be 
difficult to achieve. It is important to isolate each effect on 



the gravitational waveform in order to find general trends. If 
independent information about the source (a gravitational- 
wave burst, or optical, neutrino or gamma-ray signals in the 
case of core collapse; a gravitational-wave chirp in the case 
of a binary merger) can distinguish between isolated neutron 
star formation and binary merger, then this could be used to 
constrain the interesting range of several parameters. Thus, 
a combined filter that includes both a pre-formation char- 
acteristic signal and several damped pulsation modes will 
enhance the total signal to noise ratio (see Kokkotas, Apos- 
tolatos & Andersson, 2001, for extraction of the character- 
istic parameters of a damped monochromatic gravitational 
waves, through matched filtering). 

In the core collapse simulations by Dimmelmeier et al. 
(2001, 2002a, 2002b) the quasi-periodic gravitational waves 
emitted during core collapse were found to have frequencies 
less than roughly 1.1 kHz and our computed frequencies of 
I — and I = 2 modes agree with this, as discussed in 
Section 14.11 Such frequencies are still within the range of 
current laser-interferometric detectors. However, the above 
frequencies are typical only for certain EOSs. The funda- 
mental 1 = 2 /-mode frequency of the nonrotating model 
of sequence A is at the lower end of the corresponding fre- 
quency range for IAMq models, constructed with a large 
sample of different realistic EOSs, which range from ~ 1.35 
kHz (for extremely stiff EOSs) to ~ 3.6 kHz (for extremely 
soft EOSs), see e.g. Andersson & Kokkotas (1998). A similar- 
range of frequencies, for different EOSs, exists for the fun- 
damental I = mode, for stars of IAMq mass. Based on the 
empirical relations constructed in Section 14.11 we estimate 
that, depending on the stiffness of the high-density EOS, 
the frequency range of quasi-periodic gravitational waves 
during core collapse could be as low as 0.65 to 1.35 kHz 
for extremely stiff EOSs and as high as 1.8 to 3.6 kHz for 
extremely soft EOSs. The higher frequency range would be 
accessible for gravitational-wave detection only with detec- 
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Figure 12. Profiles of the internal energy e showing the propa- 
gation of shocks, generated due to radial pulsations at the surface 
of the star. At t = 0.14ms, a first shock is shown to propagate 
away from the star. At t = 0.92ms, a second shock is propagating 
in the rarefaction region created by the first shock. A shock is 
produced after each period of the fundamental quasi-radial mode 
F. Pulsational energy is carried away by these shocks and a high- 
entropy envelope forms in the equatorial region surrounding the 
star. 
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Figure 13. Damping of nonlinear pulsations for a star rotating 
at the mass-shedding limit. For a slowly rotating star the pulsa- 
tions are not damped (except by numerical viscosity). The time- 
evolution of the density is shown at half-radius in the star (for 
a better comparison, the lower curve was displaced by a fixed 
amount on the vertical axis). For an initial mode amplitude of 
1%, the nonlinear damping occurs on a dynamical timescale. 



tors such as the proposed wide-band dual sphere (Cerdonio 
et al. 2001). 

In the case of a neutron-star binary merger, the high- 
frequency quasi-periodic oscillations excited in the hyper- 
massive neutron star could last for a large number of os- 
cillation periods, since the damping due to a low-density 
envelope should be much weaker than in the core collapse 
case. Thus, the quasi-periodic signal could be enhanced sig- 
nificantly by matched filtering. 



8 DISCUSSION 

Using an axisymmetric general relativistic hydrodynamics 
code we have studied nonlinear pulsations of uniformly and 
differentially rotating neutron stars. We have performed 
time-dependent numerical simulations of a large sample of 
initial models which were slightly perturbed away from hy- 
drostatic equilibrium. The time-evolutions have been ana- 
lyzed using Fourier transforms at several points inside the 
stars, which enables the extraction of the two-dimensional 
eigenfunction for each mode. Our attention has been focused 
on two different sequences of uniformly and differentially ro- 
tating stars, for which we have obtained the pulsation fre- 
quencies for the two lowest-order quasi-radial (/ = 0) and 
quadrupole (I = 2) modes. We have found that differentially 
rotating models can reach significantly lower pulsation fre- 
quencies than uniformly rotating models of the same rest- 
mass. In addition, our simulations show that the fundamen- 
tal quasi-radial mode is split into two different sequences. 
This new mode has been most clearly observed in the fastest 



differentially rotating models. The centrifugal forces and the 
degree of differential rotation have significant effects on the 
mode-eigenfunctions. Most notably, axisymmetric nonradial 
modes acquire a nonzero density variation at the center of 
the star, due to coupling to lower-order terms in the eigen- 
function. However, both the splitting of the fundamental 
mode and the nonzero central density variations of nonra- 
dial modes could be due to the use of the Cowling approxi- 
mation. 

We have found that near the mass-shedding limit the 
pulsations are damped due to shocks forming at the surface 
of the star, when matter is shed in the equatorial region dur- 
ing each pulsation cycle. The damping rate should depend on 
the amplitude of the pulsations and on the rotation rate of 
the star. This new damping mechanism may set a small sat- 
uration amplitude for modes that are unstable to the emis- 
sion of gravitational radiation. The nonlinear development of 
pulsations when subject to this strong damping depends on 
the adopted EOS during the time-evolution. If one restricts 
the perfect fluid to remain isentropic, by assuming that the 
usual polytropic EOS holds throughout the evolution, then 
this restriction does not allow for real shocks to form and 
propagate away from the star, even though discontinuities 
appear in the fluid variables in this case as well (actually, 
such discontinuities are larger than the corresponding dis- 
continuities in the non-isentropic case, for the same mode- 
amplitude, due to kinetic energy conservation). Instead of 
propagating shocks, one obtains a behaviour that is remi- 
niscent of the "wave-breaking" of large- amplitude r-modes 
in rapidly rotating stars, in nonlinear simulations by Lind- 
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blom, Tohline & Valisneri (2001) 4 . In contrast, the correct 
approach is to use a nonisentropic EOS during the evolution, 
such as the ideal fluid EOS used in our simulations, which 
allows for physically realistic shocks to form and propagate. 
We note that in the simulations by Lindblom et al. (2001) 
the growth in the amplitude of the unstable r-mode was ac- 
celerated by a large factor, leading to the "wave-breaking" 
and the sudden destruction of the mode. We conjecture that 
if one would use a nonisentropic EOS and would not accel- 
erate the natural r-mode growth time, then a balance be- 
tween the growth rate due to gravitational-radiation reac- 
tion and the damping rate due to shock-dissipation caused 
by mass-shedding could be obtained at some amplitude less 
than order unity. This discussion is relevant, of course, only 
if in very rapidly rotating stars the mass- shedding- induced 
damping limits the nonlinear amplitude of unstable modes 
before other nonlinear saturation mechanisms can dominate. 
Examples of other nonlinear saturation mechanisms are sat- 
uration due to the interaction of differential rotation with 
a magnetic field (Rezzolla, Lamb & Shapiro, 2000; Rezzolla 
et al. 2001a, 2001b), hydrodynamical instabilities of large- 
amplitude nonlinear oscillations (Gressman et al. 2002) and 
nonlinear couplings to other damped modes (Schenk et al. 
2002; Morsink 2002; Arras et al. 2003). At present, it is still 
unclear what the relation, if any, between the latter two 
mechanisms is (in the case of r-mode oscillations). Which 
nonlinear saturation mechanism sets the maximum ampli- 
tude of unstable modes could depend on the rotation rate 
and the magnetic field strength of the star. 
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